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ABSTRACT 

We measure the angle between the planetary orbit and the stellar rotation axis in the transiting 
planetary system CoRoT-1, with new HIRES/Keck and FORS/VLT high-accuracy photome- 
try. The data indicate a highly tilted system, with a projected spin-orbit angle A - 11 ± 11°. 
Systematic uncertainties in the radial velocity data could cause the actual errors to be larger by 
an unknown amount, and this result needs to be confirmed with further high-accuracy spectro- 
scopic transit measurements. Spin-orbit alignment has now been measured in a dozen extra- 
solar planetary systems, and several show strong misalignment. The first three misaligned 
planets were all much more massive than Jupiter and followed eccentric orbits. CoRoT-1, 
however, is a jovian-mass close-in planet on a circular orbit. If its strong misalignment is 
confirmed, it would break this pattern. The high occurence of misaligned systems for several 
types of planets and orbits favours planet-planet scattering as a mechanism to bring gas giants 
on very close orbits. 

Key words: Planetary systems - Techniques: radial velocities - Techniques: photometric - 
Stars: individual: CoRoT-Exo-1, CoRoT-1 



1 INTRODUCTION 

The projection on the plane of the sky of the angle between a plan- 
etary orbit and the rotation axis of its host star can be measured in 
transiting systems, using the Rossiter-McLaughlin ("RM") effect 
( |Rossiter|[T924{ |McLaughIin|[T924l l. The passage of the planet in 
front of the star produces an anomaly in the radial velocity curve 
that depends on this angle. The distribution of spin-orbit angles 
is an important clue to the formation and evolution mechanisms 
of planetary systems. Close-in gas giant planets are thought to be 
formed in a disc and migrate inwards by interaction with the disc 
( |Lin etal. 1 1996) 1, a scenario expected to produce orbits with gener- 
ally aligned spin and orbit. 

The projected spin-orbit angle has now been measured pre- 



cisely for more than a dozen extra-solar planetary systems. Start- 
ing with XO-3 ( jHebrard et al.|2008] l, a substantial fraction of these 
systems exhibit large spin-orbit misalignments l |Pont et al.||2009t 
[Johnson et al.|2009) . The most extreme case to date is the retro- 
grade orbit of H AT-P-7b (|Narita et al.|2009l|Winn et al.|2009^ . 

CoRot-llp] l |Barge et al.||2008| is the first transiting planet 
identified by the CoRoT space mission. It is an extreme represen- 
tative of the "hot Jupiter" family of extra-solar planets, with an or- 
bital period of only 1.5 days. Its mass is similar to that of Jupiter 
but, like other gas giants on very close orbits, it has a much larger 
radius {R ~ 1.4 Kj). At V = 13.6 mag, the parent star lies at the 
faint end of the magnitude range for presently known planet hosts, 
so that follow-up observations require large telescopes. As part of 



* Based on observations obtained at the W.M. Keck Observatory and the 
European Southern Observatory 



Originally known as CoRoT-Exo- 1 b. The naming convention for CoRoT 
planets has been subsequently modified. 
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Figure 1. Photometric data from FORS/VLT in the B (squares) and R (cir- 
cles) filters, with the best-fit transit models. 
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Figure 2. Radial velocity data around the phase of transit. The solid line 
shows the best-fitting models. The dashed line shows the best fit with i = 
and Vrot sin/j = 5.2 kms"' imposed. 



Date [HJD] 


Relative flux 


uncertainty 


Filter 


2454524.556083 


0.9996 


0.0004 


R 



Table 2. Photometric data from FORS in B and R (full table available elec- 
tronically). 



Date [HJD] 


RV [ms-'] 


a-Rv [ms '] 


2454839.930643 


-2448.81 


17.94 



Table 1. Radial velocity measurements from HIRES/Keck and HARPS (full 
table available electronically). 



the ground-based follow-up of the CoRoT planet search, new pho- 
tometric and spectroscopic observations have been gathered on this 
target, with the objective of refining the system parameters, no- 
tably the planet size, of detecting any further longer-period planet 
in the system, and of measuring the spin-orbit angle through the 
RM anomaly. We have monitored one transit with the FORS cam- 
era on the VLT in two passbands and measured the spectroscopic 
transit with the HIRES spectrograph on the 10-m Keck telescope. 
An analysis of the photometric transit curve to constrain the planet 
size was already presented by |Gillon et al.| l |2009^ . In this study 
we concentrate on the determination of the spin-orbit angle from 
the HIRES data, using the FORS data to contrain the other transit 
parameters. 



2 OBSERVATIONS 
2.1 Photometry 

The transit was measured in high-accuracy photometry with the 
FORS camera on the VLT (ESO, Chile) in the R and B filters on 
February 27, 2008. The observing strategy and technical details are 
identical to those described in Pont et al. (2007 1. The reduction was 
described in |Gillon e t al. (2009Jl, and the data is given in Table |2] 
We chose to alternate between the R and B filter because the differ- 
ence in transit lightcurve shape between red and blue wavelengths, 
caused by the wavelength dependence of the amount of stellar limb 
darkening, offers a useful complementary constraint on the transit 
impact parameter, in complement to the shape of the transit ingress 
and egress. The mean time interval between successive measure- 
ments is 78 seconds. The dispersion of the residuals is 0.56 mmag 
in R and 0.52 mmag in B, close to the photon-noise limit. 



2.2 Radial velocity 

The radial velocity data comes from the initial follow-up spec- 
troscopy of CoRoT planet candidates with the SOPHIE and 
HARPS spectrographs, presented in [Barge et al.| ( |2008[ l, and a 
measurement of the spectroscopic transit with the Keck 1 tele- 
scope. In addition, we used radial velocity measurements obtained 
in 2008 with HARPS (ESO Prog. 082C-0.312) to constrain the 
radial-velocity orbit outside the time of transit. 

On 2009 January 8, we used Keck I and its HIRES spectrom- 
eter to observe the RM effect of CoRoT- 1 . We used the same instru- 
ment setup to measure precise RVs as the California Planet Search 
group (G. Marcy. priv. comm.). An iodine gas absorption cell is 
placed in the optical path to provide a precise wavelength calibra- 
tion. We obtained 13 spectra during the transit with an exposure 
time of 900 seconds, a slit width of 0.86 arcsec yielding a resolving 
power of R~45000, and a typical S/N ratio of « 40 : 1 in the iodine 
region. Precise differential RVs were computed with our Austral 
iodine code ( |EndI et al.|2000| l. Table [T] gives the HIRES/Keck and 
new HARPS radial velocity data. 



3 ANALYSIS 

The procedure to infer physical parameters from observations of 
transiting planetary systems has now become standard, and de- 
tailed descriptions can be found in the recent literature on the topic 
(see Winn 2009 for a recent review). A Monte-Carlo Markov 
Chain (MCMC) algorithm is used to integrate the posterior proba- 
bility distribution for the model parameters given the observed data 
and other constraints. Correlated systematics in the data are ac- 
counted for by adapting the merit function to include the effect of 
red noise. Our implementation is described in more details in Pont 



et al. 1 2009 1 and only briefly summarized here. We use the ,Mandel 
& AgoIlpOO^) and jOhta et"ar]j2005^ formahsms for the transit and 



RM effect model curves. The limb-darkening parameters are fixed 
to the values in |Claret| ( [2000^ . Stellar evolution models iGirardT] 
|et al.|2000) are used to constrain the host star's mass and radius. 
We add a factor in the merit function to take the correlation of the 
noise into account, using a set of cr,. parameters ( |Pont et al.|2006[ l, 
one for each instrument. We use the MetropoUs-Hasting jump re- 
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jection criterion for the MCMC chain, and adopt a chain length of 
10^. We verified that the correlation length for all fitted parameters 
is much shorter than the scale length. 

Based on the residuals around the model curves, we find that 
systematics dominate photon-noise uncertainties in the radial ve- 
locity data. As discussed in jPont et al.H2006) , it is important to take 
the systematics and their time correlation into account to obtain re- 
liable parameters and uncertainties. The effect of the systematics 
depends mainly of their correlation in time. There are not enough 
measurements in the radial-velocity transit sequences to measure 
the correlation length from the data itself. At the two extremes, we 
can either assume that the systematics are dominated by timescales 
very different from the transit duration, or that the systematics are 
dominated by timescale comparable to the transit duration. If, fol- 
lowing |Winn et al.| (j2008j, we note /? the ratio between the ef- 
fective uncertainties to be used in the fit to the photon-noise un- 
certainties, in the first case we have fi ^ cTrcs/Cpn, in the second 
j8 ^ n''-cr,.es/iTpn> where cTpn is the photon-noise uncertainty on sin- 
gle data points, a^^^ the r.m.s. of the residuals, and n the number 
of measurements during the transit. For our standard solution, we 
use the first value (negligible time correlation of the radial velocity 
systematics). In our experience, systematic errors for radial velocity 
measurements are less strongly correlated in time than photometric 
errors. We also examine the effect of a higher correlation. 

We fit for the following system parameters: four orbital pa- 
rameters: Tu (epoch of transit), Vo (centre-of-mass velocity), K 
(velocity semi-amplitude), ; (orbital inclination); mass and radius 
of the star Ms,Rs\ mass and radius of the planet Mp,Rp; and pro- 
jected spin-orbit angle A. We set the orbital period to the value 
found from CoroT photometry, P = 1.5089557 days, and the or- 
bital eccentricity to zero. The best constraint on eccentricity comes 
from the timing of the secondary eclipse, which indicates a negli- 
gible eccentricity (e cos o) < O.OI, |GiIlon et al.|2009| l. We use a flat 
prior in age and metallicity for CoRoT-1, as well as the tempera- 
ture and stellar rotation velocity measured in Barge et al.,p008| : 
Teff = 5950 ± 50 K, V,„t sin/, = 5.2 ± I.O km s"'. 



4 RESULTS AND DISCUSSION 
4.1 System parameters 

The best values and central 68% confidence intervals for the pa- 
rameters of the CoRoT-I system given by the MCMC integration, 
in the case assuming uncorrelated radial-velocity systematics, are 
the following: orbital elements K =198 ± 9 ms"', / = 86.7 ± 0.6°, 
a = 0.0261± 0.0005 AU, T,, = 2454524.6231± 0.0002 BJD, im- 
pact parameter b = 0.35± 0.08, star's mass and radius Mj = 1.03 ± 
0.06 Mo, i?., =1.14 ±0.03 Ro, planet's mass and radius Mp = I.I3± 
0.07 Mj, Rp = 1.48 ± 0.06Rj. We find i = 77 ±11° for the pro- 
jected spin-orbit angle. Figures[T]and|2]compare the data and best- 
fit model. Figure |2] also shows, for comparison, the best-fit model 
of the spectroscopic transit with an aligned orbit. 

The results on the mass and radius of the star and planet are 
comparable to previous determinations. Since we use the same val- 
ues of the spectroscopic parameters of the host star as [Barge et al.| 
( [2008^ , a significant update will depend on an improved determi- 
nation of the spectroscopic parameters. 

We discuss below the contraint on the spin-orbit angle of the 
system. 
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Figure 3. Probability distribution of the projected spin-orbit angle accord- 
ing to the MCMC integration. 

4.2 Spin-orbit misalignment 

Figures [3] shows the posterior probability distribution function for 
the projected spin-orbit angle. Our MCMC integration on all avail- 
able data and constraints favour an angle peaking at A = +11° 
with a probability distribution well described by a Gaussian with 
(T ^ 11°. In this solution, the star is tilted sideways relative to the 
planetary orbit, with the planet crossing only the receding side of 
the star. 

Two factors have an important influence on the values derived 
for the spin-orbit angle: the measured stellar rotation velocity, and 
the radial velocity measurement systematics during the transit. 

Rotation velocity of CoRoT-1 

A rotation velocity of Vjo, sin /, = 5.2± 1 .0 km s"' is measured from 
the broadening of the HARPS spectral cross-correlation function 
(where /, is the angle between the line-of-sight and the rotation axis 
of the star, the third angle specifying the spin-orbit configuration, 
together with A and planet orbital angle i). Since this constraint 
has an important influence on the output value for the spin-orbit 
angle, it is worth examining more closely. A lower value of the 
rotation would make the relatively flat velocity curve during the 
transit compatible with lower values for the spin-orbit angle. The 
measurement of stellar rotation from the broadening of the cross- 
correlation function has a long pedigree in radial- velocity studies, 
and the calibration for HARPS is well established. An analysis of 
the HIRES spectra also indicate significant rotation, Vrot sin /, 
6 ± 1 kms"'. We therefore consider the determination of the star's 
projected rotation velocity to be robust. 

If the a priori constraint on Vioi sin i is relaxed, the MCMC 
integration still converges towards highly tilted spin-orbit angles, 
with lower values for the rotation velocity in the 2-3 km s"' range. 
Therefore, the spectroscopic data, especially the Keck velocity 
curve, favour a high spin-orbit tilt even without external constraint 
on the stellar rotation velocity. 

Red noise in the radial velocity data 

As discussed in the case of HD 80606 in [Pont et al.| (2009), the 
posterior probability distribution for the spin-orbit angle depends 
on the choice of the assumptions on the amplitude and timescale of 
the time-correlated noise in radial velocity during the transit. 

If is of the same order as the expected RM effect (~ 50 
ms '), and correlated in time over similar timescales, then obvi- 
ously the constraint on the spin-orbit angle becomes unreliable. 

Since the time series is not much longer than the duration of 
the transit, it is not possible to determine the correlation timescale 
on the data itself, but the point-to-point scatter of the HIRES data 
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System 




e 


circular and aligned 






HD 189733b 


-1.4± 1.1 


<0.03 


HD 209458b 


0.1 ± 2.4 


< 0.03 


HAT-P-lb 


3.7 ±2.1 


< 0.10 


CoRoT-Exo-2b 


7.2 ±4.5 


< 0.10 


HD 149026b 


1.9 ±6.1 


< 0.18 


TrES-2b 


-9.0 ± 12.0 


< 0.08 


circular and tilled 






CoRoT-lb 


77± 11 


< 0.01 


HAT-P-7b 


182.5 ± 9.4 


~ 


eccentric and aligned 






HAT-P-2b 


1.2 ± 13.4 


0.52 


HD 17156b 


9.4 ± 9.3 


0.67 


eccentric and tilted 






WASP-14b 


-33 ± 7 


0.09 


XO-3b 


37.3 ± 3.7 


0.23 


HD 80606 


39 +2** 


0.93 



Table 3. Transiting exoplanets with significant projected spin-orbit angle 
determinations (formal en'ors smaller than 20°). Angles from the compi- 
lation in |Fabrycky & Winn (2009i, updated with Johnson et al. (2009); 
|Winn et al.H2009^ and this paper. Eccentricities from |Madhusudhan & Winn| 
\2Q0% except |Joshi et al.H2509t for WASP-14. 

seems to indicate correlated systematics on the 10-20 ms"' scale 
at most. This is similar to the scatter that we found on other targets 
from the CoRoT follow-up observed with the same instrument. 

We have repeated the MCMC integration using a correlation 
timescale comparable to the duration of the transit. Because the 
posterior probability distributions are relatively well-behaved and 
unimodal (see for instance Fig. [3] for the spin-orbit angle), the ef- 
fect is roughly equivalent to multiplying the parameter uncertain- 
ties by a factor n''^. With extreme assumptions on the correlation 
of the noise (n = 10), the result of a spin-orbit misalignment is still 
significant, but at a reduced level of ~ 2cr. 

Figure |2] compares the radial velocity data with the expected 
behaviour in case of an aligned orbit, assuming Vjotsin/, = 5.2 
kms"'. Explaining the data with such a model would require not 
only high systematics, but a remarkable amount of compensation, 
that we consider unlikely. However, the four data points outside the 
transit do now follow the model closely, which indicates that we do 
not have a tight control over the systematics at the point. 

One important argument in favour of a strong misalignment 
is the result from the HARPS spectroscopic transit (Bou chy et al.| 
p()09). Although the HARPS data is more noisy, it also shows a 
radial velocity constantly lower than the Keplerian orbit during the 
transit, indicating a planet crossing only the receding half of the 
star. 

Precise and high-cadence radial velocity measurements for an 
object fainter than 13th magnitude in the visible are difficult even 
with the best instruments and largest telescopes. As for previous 
cases of measured spin-orbit misalig nment, XO-3 jHebrard et al. 
2008;Winn etal.'2009l, HD 80606 ( Moutou et al.|2009[|Winn etaT 
2009) and WASP-14 ( Johnson et al. 20 09 1, further observations and 
an independent confirmation of our result will be useful. 

4.3 Discussion 

The projected spin-orbit angle has now been measured with ade- 
quate accuracy for more than a dozen transiting planets. The main 
result is the strikingly frequent occurrence of high spin-orbit mis- 
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Figure 4. Distribution of projected spin-orbit angles for systems in Table|3] 



alignments. At least four systems have spin-orbit tilts in excess 
of 30 degrees (XO-3, HD 80606, WASP-14, HAT-P-7, see Intro- 
duction). Evidence for misalignment may also have been found for 
WASP- 17 (Anderso n et al.|2009[ ). According to the data and analy- 
sis presented here, CoRoT-1 joins this group, with the radial veloc- 
ity data favouring an almost polar orbit. 

Table [3] displays the systems with significant spin-orbit an- 
gle measurements, in relation with orbital eccentricities. Figure |4] 
shows this list in histogram form. Note that it is important for this 
ensemble analysis to include only conclusive determinations that 
do account for correlated noise in the uncertainties. Values with un- 
derestimated uncertainties introduce an obvious bias towards non- 
alignment. 

I Johnson et al.| ( [2009] l pointed out an interesting feature of the 
first three misafigned systems (XO-3, WASP-14 and HD 80606): 
all contained a very high-mass gas giants (4, 7 and 12 Jupiter 
masses) on markedly eccentric orbits. This suggested that misalign- 
ment was not a feature of ordinary hot jupiters, but of their higher- 
mass cousins. CoRoT-1, however, would break this pattern. It has a 
mass nearly equal to Jupiter, and a circular orbit. The same is true 
for HAT-P-7 (Narita et al.p009| [WmiT et al.|2009| >. Therefore the 
correlation between planet mass, eccentricity and misalignment is 
called into question. 

The cases of CoRoT-1 and HAT-P-7 reinforce the statisti- 
cal conclusions drawn on an earlier sample by |Fabrycky & Winn| 
(2009): the sample as a whole is compatible with a bimodal dis- 
tribution of spin-orbit angles, with some systems well aligned, and 
others with little correlation between the two angles. The distribu- 
tion illustrated by Fig.|4]is compatible with the combination of two 
sets of system, one with aligned orbits, and the other with randomly 
oriented orbits. As more systems are measured, most precise statis- 
tical properties of the distribution of orbital inclinations for close-in 
giant planets will become accessible. 

An almost polar orbit of CoRoT-lb would be very sugges- 
tive of the action of an extreme dynamical event, such as a close 
encounter with another massive planet. Formation of close-in gas 
giants by migration in a gas disc would mainly result in aligned 
orbits jLin et al.|I996^ . Gravitational interactions with other bod- 
ies can perturb the spin-orbit alignment both during and after the 
planet formation stage, either through orbital angular momentum 
exchange (e.g. Kozai mechanism) or more dramatically through 
catastrophic resonances and close encounters. The first process is, 
for instance, a plausible explanation for the extremely eccentric 
system HD 80606. It is less likely to operate on CoRoT-lb, given 
its close and circular orbit and the absence of detected stellar or 
planetary perturber. The second type of process is a more likely 
explanation. 

It has long been suggested that planet-planet scattering may 
account for the presence of close-in exoplanets l |Chatterjee et al.| 
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|2008[ [Ford & Rasio|[2008] l. High spin-orbit angles would seem 
to favour this hypothesis, which has also gained renewed support 
from other lines of evidence: multiple planetary systems tend to 
be dynamically dense ( ,Raymond et al.|[2009j , the distribution of 
planetary eccentricities resembles that expected from a scattering 
process j juric & Tremaine|2008| l, and the closest orbital distances 
conform to the scattering scenario ( |Ford & Rasio|2006^ . The case 
of CoRoT- 1 illustrates the high interest of vigorously pursuing the 
effort to measure the RM anomaly in transiting exoplanets. 
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